Phase diagram of a surface superconductor in parallel magnetic field 
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Universal phase diagram of 2D surface superconductor with generic Rashba interaction in a paral- 
lel magnetic field is found. In addition to uniform BCS state we find two inhomogeneous supercon- 
ductive states, the stripe phase with A(r) oc cos(Qr) at high magnetic fields, and a new "helical" 
phase with A(r) cx exp(iQr) which intervenes between BCS state and stripe phase at intermediate 
magnetic field and temperature. We prove that the ground state for helical phase carries no current. 
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Many efforts both theoretical and experimental have 
been spent in search of exotic nonhomogeneous super- 
conducting states beginning with a pioneering papers 
by Larkin-Ovchinnikov and Fulde-Ferrel (LOFF)i*^ pre- 
dicting the stripe state in superconductors with compet- 
ing ferromagnet interaction. Nevertheless no convinc- 
ing experimental evidence for existence of such a state 
was found till now, partially due to rather narrow ex- 
istence range of the LOFF state. Recently Barzykin 
and Gor'kov'^ did find a system where such inhomoge- 
neous superconducting phase could be prominent. It 
is a two-dimensional surface superconductor with spin- 
orbital Rashba interaction. One of possible realization 
of such a system probably was reported in very inter- 
esting experiments^, where various signatures of surface 
superconductivity with Tc ~ 90if were detected in the 
insulating WO^ doped by small amount of Na. Surface 
spin-orbit superconductivity is unusual one due to the 
absence of inversion symmetry: this results in the pres- 
ence of the the spin-orbital Rashba term^ and the chiral 
subband splitting of free electron spectrum at the surface. 
In such a superconductor condensate wave-function is a 
mixture of both singlet and triplet stateaS^, therefore 
Pauli susceptibility is not vanishing'' at T — > 0; para- 
magnetic breakdown of superconductivity in a parallel 
magnetic field is shifted towards much higher field values 
due to the formation of LOFF stated. The line of transi- 
tion from normal to (any of) superconductive state Tc{h) 
was determined in^; however, the nature of the phase in- 
tervening between uniform BCS and stripe LOFF state 
was not studied. This question is important because phe- 
nomenological theory^ predicts the possibility for a new 
helical state distinct from both BCS and stripe LOFF 
states. 

In this Letter we provide a detailed phase diagram of 
a surface superconductor in a parallel magnetic field h. 
The phase diagram turns out to be universal after nor- 
malization of the temperature and the Zeeman energy 
by the critical temperature in zero magnetic field Tco for 
two models of high (I) and low (II) electron density. The 
model I assumes a normal 2D metal with Fermi energy 
being much greater than chiral splitting. The model II 
is suited for an electron gas in field effect heterostruc- 
tures where electrons fill only the bottom of one of chiral 
bands. We demonstrate the existence of a helical state 



with order parameter A (x exp(iQr) (where Q _L h) 
and Q fish/vp in a considerable part of the phase di- 
agram, which is summarized in Fig. 1 . The line CT is the 
second-order transition line separating helical state from 
the homogeneous superconductor. Below the T point 
first-order transition between BCS and inhomogeneous 
state takes place. The line ST' is the line of soft insta- 
bility of the helical state (see below). The point S in 
the Tc{h) line is special in that here the order parameter 
symmetry is U{2) instead of usual U{1). Full details of 
our theory will be presented in a separate publication^. 

Near the surface of a crystal the translational symme- 
try is reduced and the inversion symmetry is broken even 
if it is present in the bulk. As a result a transverse elec- 
tric field appears at interface and gives rise to the rela- 
tivistic spin-orbit interaction known as the Rashba term: 

Hso — a. (7 X p • n, where a > is the spin-orbit cou- 
pling constant, n is a unit vector perpendicular to the 
surface, u = (cr^, cr^, cr^) are spin Pauli matrices. Spin 
operator does not commute with the Rashba term, thus 
spin projection is not a good quantum number. On the 
other hand the chirality operator: cr^ sini^Sp — cosi^Sp, 
does commute with the Hamiltonian, where ip-p is the 
angle between the momentum of the electron and the 
X-axis. Its eigenvalue: A = ±1 is a quantum number 
of an electron state (p, A) . The Rashba term preserves 
the Kramers degeneracy with states (p, A) and (— p. A) 
belonging to the same energy. 

In this Letter we consider the simplest model of sur- 
face superconductor: Gor'kov model for two-dimensional 
metal with the Rashba term being included^, in the limit 
apF S> Tc. The Hamiltonian written in the coordinate 
representation reads 



a = I ?A+(r) ( ^Sai3 + a aap X P 
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where m is the electron mass and a, (3 are spin indices, 
P — — iV-|- ^A{r) is the momentum operator in the pres- 
ence of infinitesimal in-plane vector-potential A — A{r). 
Zeeman interaction with a uniform external magnetic 
field h parallel to the interface and in the z-direction 
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is included. The vector-potential of such a field can be 
chosen to have only z-component, therefore it decouples 
from the 2D kinetic energy term, /is is Bohr magneton 
and g is the Lande factor. Hereafter we use a notation 
H = gfiBh/2. 

The electron operator can be expanded in the basis of 
plane waves tpair) — X^Ap c^^'^Cap. The one-particle part 
of the Hamiltonian in the momentum representation: 

Ha = Y^ c+p + a [a^f} xp\-n~H ■ ct^^"] c^p (2) 



can be diagonalized by the transformation c^p — 
Vxaipjaxp with the two-component spinor r/\{p) — 

(1, zAexp(j(^p))/\/2- Eigenvalues of the Hamiltonian (O 
corresponding to the chiralities A = ±1 are 



^xip) = p^/2m - X^a^p^ - 2apyH + H^. (3) 

In the case of high electron density and respectively 
chemical potential fi ^ ma^ ( model I) both chiral- 
ity branches are filled and the equal momentum electron 
states are split by 2apF- Fermi circles with different chi- 
ralities are split with the radii = \/2rn/i -t- m^a^ima. 
Densities of states on the two Fermi circles are almost the 
same, v± = ^ {l± , and in this paper we neglect 
the difference — v_. In the case of the low electron 
density (model II): — ma^/2 < /i < 0, the electrons fill 
the bottom of only one chiral branch A = 1, i.e. a ring. 
The two Fermi circles of the ring in zero magnetic field 
have radii pf = wa -|- ly^2mij + m?a^^ where I = ±1 
is the number specifying the inner and the outer part of 
the ring and is analogous to chirality. Density of states 
on the outer and inner part of the Fermi ring (I — ±1) 
is almost the same = in the case of narrow 

ring: vp <C oi. If magnetic field is applied, the two Fermi 
circles are displaced in opposite y-directions by a mo- 
mentum Q = zLH/vf, where vp = \/2^/m + is the 
Fermi velocity in the I model, or half width of the ring 
divided by m in the II model. The pairing interaction 
(the last term in Eq.lQJ) can be factorized in the chiral 
basis: Hint = — -j X]q^^('Z)^('?)' where the pair an- 
nihilation operator A{q) — J2pX '^^^^^^x-p-^Xp+ '^^ith 

p± = p± q/2. 

To calculate thermodynamic potential f2 = — TlnZ, 
we employ imaginary-time functional integration tech- 
nique with Grassmanian electron fields a\ p, a\,p and in- 
troduce auxiliary complex field A(r, r) to decouple pair- 
ing term Hint, cii^. The resulting effective Lagrangian 
is of the form: 

L[a, a, A, A*] = ^ axp (-9r - £x{p)) axp + 
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Below we will work within the mean-field approxima- 
tion which is controlled by the small Ginzburg number 
Gi ^ Tc/Ep. It is equivalent to the saddle-point approx- 
imation for the functional integral over A and A*, thus 
we will be studying minima of the functional r2[A,A*] 
appearing after Gaussian integration over Grassmanian 
fields. 

Near the normal-superconducting phase transition the 
order parameter A(r) is small and may be expanded 
in powers of A(r) and its gradients. We consider the 
order parameter as a superposition of the finite number of 
harmonics, A(r) = ^ ■ Aq. (r) exp (iQ,;r), where Aq. (r) 
are slowly varying in 2D space functions, and derive the 
corresponding Ginzburg-Landau functional: 
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The coefficient a(Q) includes Cooper loop diagram with 
transferred momentum Q: 
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Gl{uj,p+)Gl{-uj,-p-L) (6) 



where in the I model the normal-state Green function in 
in-plane magnetic field H <ti app is {^iP) = {i^ ~ 
^ — XH shiipp)^^ , and ^ = p^/2m~ Xapp — /it is assumed 
to be small compared to app, whereas in the II model 
H <C ma^ and ^; = vp{l'K — mvp) is assumed to be 
small compared to ma^ (tt is related to the center of the 
ring). Hereafter we present results for high density (I 
model), but results for low density are identical. The 
condition ming a{Q) — determines the second-order 
transition line (if /3 > 0) between the normal metal and 
the superconductor: 
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hereafter Hx — XH + vpQ/2. Depending upon if, the 
maximum in Eq.ljJJ) is attained either at Q = or at 
nonzero ±|(5|. Tc{H) line is shown on Fig^where both 
Tc and H are normalized by the critical temperature 
at zero magnetic field Tco = 2ujp)exp{—l/i'U + G)/tt, 
where C = 0.577 is the Euler constant. Near C point 
vpQ{H) - - Hi. In the hmit H/T^o ^ oo we re- 

cover the asymptotics of Tc{H) linei, and find vpQ{H) = 

_4rp4 

acj^a ■ Note that Q = 2H/vp is the momen- 
±1 Fermi surfaces in parallel 
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tum splitting of two A 
magnetic field. Lifshitz point C separates Q — and 
Q ^ solutions on the Tc{H) fine; it is the end point of 
the second order phase transition between the two super- 
conducting phases. At _ff > H]^ inhomogeneous super- 
conductive phase is formed below Tc{H) line. No more 
than two harmonics contribute to A(r) just below Tc{H) 
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FIG. 1: Phase diagram that shows: superconducting phase 
transition line Tc{H) (bold solid) and two second order phase 
transition lines in the clean case: CT line between homoge- 
neous (BCS) and helical (h.) state and ST' line of stabil- 
ity of helical state. Short-dashed line going downwards from 
the point T marks the absolute limit of stability of the BCS 
state. Dotted line shows the physical Tbkt (H) line for values 
Tco/ep = 0.02 and a/vp = 0.34. 

line: A(y) = A+e^^^' -I- A_e~*^y. Below Tc{H) density 
of the thermodynamic potential f2 is lower in the super- 
conductive state than in the normal one by the amount: 

Qsn = «|A|2 + /5,|A|4 + /3,(|A+|2 - \A^\y, (8) 

where jAp = jA+p + jA.p. We find coefficients Ps,a 
using standard diagram expansion around normal state. 
At the symmetric point S where l3a{Tc{H),H) = the 
free energy ((HJ is invariant under U{2) rotations of the 
order parameter spinor (A_|_,A_). At H < Hs we find 
/3a < and the free energy at T < Tc{H) is minimized by 
the choice of either A+ = or A_ = 0, both correspond 
to helical state. At iJ > Hs we find /3q > and in the 
free energy minimum |A_|_| = |A_|, which corresponds 
to the LOFF-like stripe phase with A{y) cx cos{Qy). At 
lower temperatures in this phase A.(y) contains higher 
harmonics and is time reversal symmetric. In the heli- 
cal state with only one harmonics A(y) = Ae*'^^, the 
thermodynamic potential reads: 

Qhei v/(^ + *^Asin(^)2 + A2^ + ^.(9) 

In equilibrium the stationary conditions dflhei/dA ~ 
and dilhei/dQ — are satisfied and they can be found 
explicitly: 

where r{Hx,uj) = y/uj'^ + {\Hx\ + A)^, the Jacoby mod- 
ulus k — A\Hx\/r{Hx,u;) and the function f{Hx,uj) 



is defined through the Jacoby complete elliptic integrals 
of the first and the second kind, 

f{H.) = f (c.^ + Hi + A^)^|^ - K^A, -)E(fc)) . 
Hx V r{Hx,uj) J 

ill) 

We prove by direct microscopic calculation that js = 
^di}/dQ, therefore the equilibrium state carries no su- 
per current. 

Minimum of the thermodynamic potential (j^ over 
A can be expanded in series of small Q: ^hei{Q) — 
^hei{Q) + ^nfiQ"^ + bQ^ + cQ^, c > 0. The condi- 
tion = 0, & > determines the second order Lifshitz 
transition line >CT, which ends at the point T, where co- 
efficient & = changes sign. We compute the coordinates 
of T point using Eas. (|9ll()|l . At lower temperatures & < 
and first-order transition out of BCS state occurs. 

The domain of helical state local stability was deter- 
mined via consideration of the thermodynamic potential 
variations due to weak static modulation of the form 
(5A(r) = exp(-ig?/) -t- u^+2Qexp(i(q -I- 2Q)y) (the 
presence of two Fourier harmonics in the perturbation 
is due to inhomogeneity of the background helical state): 
SVtsv — v'^A{q)v, where v = (Sv-q, 5v*^2q) ^-^d 

X J__ / Gap_Ga-j9+ Fxp_Fx^p^ \ 

(12) 

The matrix A has two eigenvalues ei{q) < £2(9)- We de- 
fine the helical state metastability line ST' as a collection 
of points where one mode Sv becomes energetically favor- 
able: mingei(g) =0. We solve numerically four equations 
simultaneously: two gap equations (|10|l that determine 
equilibrium A and Q, together with the two equations 
dq^iil) = and ei{q) — 0. By means of expansion of 
the Ginzburg-Landau functional up to the terms of order 
|A-|-|®, we checked that the next-order "anisotropic" term 
in the expansion ||SJ) is of the form e(| A+p — | A_ p)^, with 
e > 0. This fact ensures that the phase transition out of 
the helical state is of the second order (at least, near the 
T,{H) line). 

We calculated electromagnetic response function 
Sja/SAf3 = —^n'f^ for helical state using standard di- 
agram methods and found that = 4^|^. Thus on 
the Lifshitz line CT there is no linear supercurrent in 
the direction perpendicular to the magnetic field. The 
component nf ^ does not vanish anywhere in helical state 
region and is of the order of of the BCS state. This 
is in contrast with the classical LOFF problem, where 
Us was shown to vanish in the whole helical state; the 
difference is probably due to the fact that in our problem 
the direction of Q is fixed by an applied field h, while 
for the case of ferromagnetic superconductor it is arbi- 
trary. The obtained behaviour of n"^ tensor indicates 
strongly anisotropic electromagnetic response of surface 
superconductor near the Lifshitz line CT. 

So far we discussed the clean case; below we demon- 
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strate that sufRciently high concentration of nonmagnetic 
impurities suppress both hehcal and stripe states. Con- 
sider impurities with weak short-range potential, char- 
acterized by to elastic scattering time t. By means of 
diagram techniqusiS, we calculate the coefficient a{Q) 
in the Ginzburg-Landau expansion ^ which is the 
electron-electron vertex in the Cooper channel in the 
presence of nonmagnetic impurities. It is given by a 
sum of ladder diagrams which are an alternating se- 
quence of blocks of two normal metal Green functions 
Gx = {iu! ~ S, — XHsinipp + ^sgnw)"^ and an impu- 
rity line. In every block momenta on the upper and 
lower lines are opposite whereas the chiralities are the 
same. Tc{H) line is found from ming a{Q) = 0, where 

a{Q) = 1/U — TrvT maxg J2lj>o ^ ('^j Q) ' with the 
Cooper kernel K given by 



(13) 



where I^^ 



= II ± II (7 = 0,1,2) are func- 
tions of {T,H,T,Q): = {Cb^ +Hiy^^^ /4t, !{ = 

I (cl^/o - 4^) /iJA, II = *c^/l/-ffA, Co ^ u: + ^. At 
H = time- reversal symmetry is recovered and Ea. H13f) 
simplifies to K = 2/uj independently on disorder, in 
agreement with Anderson theorem. We evaluate numer- 

ically a — d^a{Q)/dQ along the the transition line 
a{T,H,Q) = 0. At rTco/ft < 0.11 we found a" > 
at any H, i.e. both stripe and helical state disappear 
from the phase diagram at t < h/9Tco. Similar condi- 
tion for usual LOFF state is more stringent, ciM>. In 
the dirty limit tTcq ^ 1 the kernel K simplifies to 
K = 2/{oj + 2H^T + vIQ^t/4). From this form of K 
one easily concludes that the paramagnetic critical field 
grows with increase of disorder as Hp = 7rTco/4re'^, cf. 
similar result in^"*. 

The phase diagram Fig^was obtained neglecting small 
parity-breaking term of the order a/vp <C 1 in the 
thermodynamic potential Sfl — rjQ with 77 = 

—avT^^f(H,uj), where function / is defined for the 
clean case in Eq. ()ll|l . Taking this term into account while 



minimizing 11, one finds that uniform BCS state trans- 
forms into a weakly helical state (predicted phenomeno- 
logically in^) with small wave-vector Q w and with- 

out supercurrent (homogeneous state would carry super- 
current, as was found inii, but it is not the ground-state). 
The line of 2nd order transition CT broadens then into 
a sharp cross-over region between two helical states with 
small and large values of Q. Within this crossover region, 
the superfluid density tensor is strongly anisotropic, with 
/ri^^ ^ {a/vpY^'^ ■ In the dirty limit long-wavelength 
helical modulation is present everywhere in supercon- 
ducting state; near the transition line its wave-vector 

We have calculated Tc{H) line within mean-field ap- 
proximation (MFA), those accuracy is usually of or- 
der Tc/ep for clean 2D superconductor: actual transi- 
tion is of Berezinsky-Kosterlitz-Thouless vortex depair- 
ing type, and is shifted downwards in temperature by 
about T^/ep. In our system fluctuations are enhanced 
strongly around C and S points. Near the C point it is 
due to smallness of n|^, and the enhancement factor is 
of the order of a/ n^^ /nV ~ (a/wj?)"^/'^. In the vicinity 
of the point S fluctuations are enhanced due to extended 
U{2) symmetry of the order parameter; 2D renormaliza- 
tion group calculation shows that U{2) fluctuation modes 
shift actual by AT^ ~ ^T'^/ep) \ogl3s/(3a downwards 
at f3a ^ Ps- In result, phase transition line Tc{H) is de- 
formed in the vicinities of C and S points, as shown on 
Fig.l. 

In conclusions, we have demonstrated an existence of 
three different superconductive states (two helical states 
and stripe state) in clean surface superconductor in par- 
allel magnetic fleld, and have located transition lines be- 
tween them. Strong disorder eliminates short-wavelength 
helical and stripe states, whereas long-wavelength heli- 
cal state survives. We thank P. M. Ostrovsky and M. 
A. Skvortsov for many useful discussions. This research 
was supported by SCOPES grant, RFBR grant 01-02- 
17759, Russian ministry of science and Program "Quan- 
tum macrophysics" of RAS. O. V. Dimitrova is grateful 
to Dynasty Foundation for financial support. 
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